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ABSTRACT 

Pupil mapping is a promising and unconventional new method for high con- 
trast imaging being considered for terrestrial exoplanet searches. It employs two 
(or more) specially designed aspheric mirrors to create a high-contrast amplitude 
; profile across the telescope pupil that does not appreciably attenuate amplitude. 

As such, it reaps significant benefits in light collecting efficiency and inner work- 
ing angle, both critical parameters for terrestrial planet detection. While much 
has been published on various aspects of pupil mapping systems, the problem of 
sensitivity to wavefront aberrations remains an open question. In this paper, we 
present an efficient method for computing the sensitivity of a pupil mapped sys- 
tem to Zernike aberrations. We then use this method to study the sensitivity of 
a particular pupil mapping system and compare it to the concentric-ring shaped 
pupil coronagraph. In particular, we quantify how contrast and inner working 
angle degrade with increasing Zernike order and rms amplitude. These results 
have obvious ramifications for the stability requirements and overall design of a 
planet-finding observatory. 



Subject headings: Extrasolar planets, coronagraphy, Fresnel propagation, diffrac- 
tion analysis, point spread function, pupil mapping, apodization, PIAA 
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1. Introduction 

The impressive discoveries of large extrasolar planets over the past decade have inspired 
widespread interest in finding and directly imaging Earth-like planets in the habitable zones 
of nearby stars. In fact, NASA has plans to launch two space telescopes to accomplish this, 
the Terrestrial Planet Finder Coronagraph (TPF-C) and the Terrestrial Planet Finder In- 
terferometer (TPF-I), while the European Space Agency is planning a similar multi-satellite 
mission called Darwin. These missions are currently in the concept study phase. In addition, 
numerous ground-based searches are proceeding using both coronagraphic and interferomet- 
ric approaches. 

Direct imaging of Earth-like extrasolar planets in the habitable zones of Sun-like stars 
poses an extremely challenging problem in high-contrast imaging. Such a star will shine 10 10 
times more brightly than the planet. And, if we assume that the star-planet system is 10 
parsecs from us, the maximum separation between the star and the planet will be roughly 
0.1 arcseconds. 

Design Concepts for TPF-C. For TPF-C, for example, the current baseline design 
involves a traditional Lyot coronagraph consisting of a modern 8th-order occulting mask 
(see, e.g., Kuchner et al. (2004)) attached to the back end of a Ritchey- Chretien telescope 
having an 8m by 3.5m elliptical primary mirror. Alternative innovative back-end designs 
still being considered include shaped pupils (see, e.g., Kasdin et al. (2003) and Vanderbei 
et al. (2004)), a visible nuller (see, e.g., Shao et al. (2004)) and pupil mapping (see, e.g., 
Guyon (2003) where this technique is called phase-induced amplitude apodization or PIAA). 
By pupil mapping we mean a system of two lenses, or mirrors, that takes a flat input field 
at the entrance pupil and produce an output field that is amplitude modified but still flat in 
phase (at least for on- axis sources). 

The Pupil Mapping Concept. The pupil mapping concept has received considerable 
attention recently because of its high throughput and small effective inner working angle 
(IWA). These benefits could potentially permit more observations over the mission lifetime, 
or conversely, a smaller and cheaper overall telescope. As a result, there have been numerous 
studies over the past few years to examine the performance of pupil mapping systems. In 
particular, Guyon (2003); Traub and Vanderbei (2003); Vanderbei and Traub (2005); Guyon 
et al. (2005) derived expressions for the optical surfaces using ray optics. However, these 
analyses made no attempt to provide a complete diffraction through a pupil mapping system. 
More recently, Vanderbei (2006) provided a detailed diffraction analysis. Unfortunately, this 
analysis showed that a pupil mapping system, in its simplest and most elegant form, cannot 
achieve the required 10^ 10 contrast; the diffraction effects from the pupil mapping systems 
themselves are so detrimental that contrast is limited to 10~ 5 . In Guyon et al. (2005) and 
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Pluzhnik et al. (2006), a hybrid pupil mapping system was proposed that combines the pupil 
mapping mirrors with a modest apodization of oversized entrance and exit pupils. This 
combination does indeed achieve the needed high-contrast point spread function (PSF). In 
this paper, we call such systems apodized pupil mapping systems. 

A second problem that must be addressed is the fact that a simple two-mirror (or 
two-lens) pupil mapping system introduces non-constant angular magnification for off-axis 
sources (such as a planet). In fact, the off-axis magnification for light passing through a small 
area of the exit pupil is directly proportional to the amplitude amplification in that small 
area. For systems in which the exit amplitude amplification is constant, the magnification 
is also constant. But, for high-contrast imaging, we are interested in amplitude profiles that 
are far from constant. Hence, off-axis sources do not form images in a formal sense (the 
"images" are very distorted.) Guyon (2003) proposed an elegant solution to this problem. 
He suggested using this system merely as a mechanism for concentrating (on-axis) starlight 
in an image plane. He then proposed that an occulter be placed in the image plane to remove 
the starlight. All other light, such as the distorted off-axis planet light, would be allowed to 
pass through the image plane. On the back side would be a second, identical pupil mapping 
system (with the apodizers removed), that would "umap" the off-axis beam and thus remove 
the distortions introduced by the first system (except for some beam walk — see Vanderbei 
and Traub (2005)). 

Sensitivity Analysis. What remains to be answered is how apodized pupil mapping 
behaves in the presence of optical aberrations. It is essential that contrast be maintained 
during an observation, which might take hours during which the wavefront will undoubt- 
edly suffer aberration due to the small dynamic perturbations of the primary mirror. An 
understanding of this sensitivity is critical to the design of TPF-C or any other observatory. 
In Green et al. (2004), a detailed sensitivity analysis is given for shaped pupils and vari- 
ous Lyot coronagraphs (including the 8th-order image plane mask introduced in Kuchner 
et al. (2004)). Both of these design approaches achieve the needed sensitivity for a realizable 
mission. So far, however, no comparable study has been done for apodized pupil mapping. 
One obstacle to such a study is the considerable computing power required to do a full 2-D 
diffraction simulation. 

Aberrations Given by Zernike Polynomials. In this paper, we present an efficient 
method for computing the effects of wavefront aberrations on apodized pupil mapping. We 
begin with a brief review of the design of apodized pupil mapping systems in Section 2. We 
then present in Section 3 a semi-analytical approach to computing the PSF of systems such 
as pupil-mapping and concentric rings in the presence of aberrations represented by Zernike 
polynomials. For such aberrations, it is possible to integrate analytically the integral over 
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azimuthal angle, thereby reducing the computational problem from a double integral to a 
single one, eliminating the need for massive computing power. 

In Section 4, we present the sensitivity results for an apodized pupil mapping system 
and a concentric ring shaped pupil coronagraph, and compare the results. 



In this section, we review the apodized pupil mapping approach and introduce the 
specific system that we study in subsequent sections. It should be noted that this apodized 
pupil mapping design may not be the best possible. Rather, it is merely an example of such 
a system that achieves high contrast. Other examples can be found in the recent paper 
by Pluzhnik et al. (2006). Our aim in this paper is not to identify the best such system. 
Instead, our aim is to develop tools for carrying a full diffraction analysis of any apodized 
pupil mapping system. 



We begin by summarizing the ray-optics description of pure pupil mapping. An on-axis 
ray entering the first pupil at radius r from the center is to be mapped to radius f = R(r) 
at the exit pupil (see Figure 1). Optical elements at the two pupils ensure that the exit ray 



or, sometimes, negative and decreasing. In either case, the function has an inverse that 
allows us to recapture r as a function of f: r = R(r). The purpose of pupil mapping is to 
create nontrivial amplitude profiles. An amplitude profile function A{f) specifies the ratio 
between the output amplitude at f to the input amplitude at r (in a pure pupil-mapping 
system the input amplitude is constant). Vanderbei and Traub (2005) showed that for any 
desired amplitude profile A(r) there is a pupil mapping function R{r) that achieves it (in a 
ray-optics sense). Specifically, the pupil mapping is given by 



2. Review of Pupil Mapping and Apodization 



2.1. Pupil Mapping via Ray Optics 



is parallel to the entering ray. The function R(r) is assumed to be positive and increasing 




(1) 



Furthermore, if we consider the case of a pair of lenses that are planar on their outward-facing 
surfaces, then the inward-facing surface profiles, h(r) and h(r), that are required to obtain 
the desired pupil mapping are given by the solutions to the following ordinary differential 
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equations: 



and 



dh r — R(r) 



dr 



r) = W = =, HO) = z, (2) 



\n 



l\^ 2 + ^(r-R(r)) 2 



§(0 = , m ~ f *(o) = o. (3) 

Here, n ^ 1 is the refractive index and z is the distance separating the centers (r = 0, f = 0) 
of the two lenses. 

Let S(r, r) denote the distance between a point on the first lens surface r units from the 
center and the corresponding point on the second lens surface f units from its center. Up to 
an additive constant, the optical path length of a ray that exits at radius f after entering at 
radius r = R(r) is given by 

Q (f) = S(R(f),f) + \n\(h(f) - h(R(f))). (4) 

Vanderbei and Traub (2005) showed that, for an on-axis source, Qo(r) is constant and equal 
to -(n - Did. 1 



2.2. High-Contrast Amplitude Profiles 

If we assume that a collimated beam with amplitude profile A(f) such as one obtains 
as the output of a pupil mapping system is passed into an ideal imaging system with focal 
length /, the electric field E(p) at the image plane is given by the Fourier transform of A(f): 

oo oo 

E ^ ,r]) = ^f Cm ^ I I e ' 2m ^ A (V^+¥)dydx. (5) 

— oo — oo 

Here, Eq is the input amplitude which, unless otherwise noted, we take to be unity. Since the 
optics are azimuthally symmetric, it is convenient to use polar coordinates. The amplitude 
profile A is a function of f = ^/x 2 + y 2 and the image-plane electric field depends only on 



For a pair of mirrors, put n = — 1. In that case, z < as the first mirror is "below" the second. 
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image-plane radius p = a/£ 2 + i] 2 : 



OO 2-7T 



E{p) 



Xif 



1 



cos 



V-*)A(r)rd0dr 



(6) 



o o 



oo 




-27T-y A(f)fdf. 



(7) 







The point-spread function (PSF) is the square of the electric field: 



Psf(p) = |£(p)| 2 . 



(8) 



For the purpose of terrestrial planet finding, it is important to construct an amplitude profile 
for which the PSF at small nonzero angles is ten orders of magnitude reduced from its value 
at zero. A paper by Vanderbei et al. (2003a) explains how these functions are computed as 
solutions to certain optimization problems. The high-contrast amplitude profile used in the 
rest of this paper is shown in Figure 2. 



Vanderbei (2006) showed that pure pupil mapping systems designed for contrast of 10~ 10 
actually achieve much less than this due to harmful diffraction effects that are not captured 
by the simple ray tracing analysis outlined in the previous section. For most systems of 
practical real- world interest (i.e., systems with apertures of a few inches and designed for 
visible light), contrast is limited to about 10 -5 . Vanderbei (2006) considered certain hybrid 
designs that improve on this level of performance but none of the hybrid designs presented 
there completely overcame this diffraction-induced contrast degradation. 

In this section, we describe an apodized pupil mapping system that is somewhat more 
complicated than the designs presented in Vanderbei (2006). This hybrid design, based on 
ideas proposed by Olivier Guyon and Eugene Pluzhnik (see Pluzhnik et al. (2006)), involves 
three additional components. They are 

1. a preapodizer A to soften the edge of the first lens/mirror so as to minimize diffraction 
effects caused by hard edges, 

2. a postapodizer to smooth out low spatial frequency ripples produced by diffraction 
effects induced by the pupil mapping system itself, and 



2.3. Apodized Pupil Mapping Systems 
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3. a backend phase shifter to smooth out low spatial frequency ripples in phase. 

Note that the backend phase shifter can be built into the second lens/mirror. There are 
several choices for the preapodizer. For this paper, we choose Eqs. (3) and (4) in Pluzhnik 
et al. (2006) for our pre-apodizer: 



A(r)+f3A max ' 

where A max denotes the maximum value of A(r) and f3 is a scalar parameter, which we take 
to be 0.1. It is easy to see that 

• A{r)/A max < A (r) < 1 for all r, 

• A (r) approaches 1 as i(r) approaches A max , and 

• A (r) approaches as A(r) approaches 0. 

Incorporating a post-apodizer introduces a degree of freedom that is lacking in a pure 
pupil mapping system. Namely, it is possible to design the pupil mapping system based on 
an arbitrary amplitude profile and then convert this profile to a high-contrast profile via an 
appropriate choice of backend apodizer. We have found that a simple Gaussian amplitude 
profile that approximately matches a high-contrast profile works very well. Specifically, we 
used 

^4pupmap (0 = 3.35e (~/~) ^ 
where a denotes the radius of the second lens/mirror. 

The backend apodization is computed by taking the actual output amplitude profile 
as computed by a careful diffraction analysis, smoothing it by convolution with a Gaussian 
distribution, and then apodizing according to the ratio of the desired high-contrast amplitude 
profile A(f) divided by the smoothed output profile. Of course, since a true apodization can 
never intensify a beam, this ratio must be further scaled down so that it is nowhere greater 
than unity. The Gaussian convolution kernel we used has mean zero and standard deviation 

s/vioo,ooo. 

The backend phase modification is computed by a similar smoothing operation applied 
to the output phase profile. Of course, the smoothed output phase profile (measured in 
radians) must be converted to a surface profile (having units of length). This conversion 
requires us to assume a certain specific wavelength. As a consequence, the resulting design 
is correct only at one wavelength. The ability of the system to achieve high contrast degrades 
as one moves away from the design wavelength. 
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2.4. Star Occulter and Reversed System 

It is important to note that the PSFs in Figure 2 correspond to a bright on-axis source 
(i.e., a star). Off-axis sources, such as faint planets, undergo two effects in a pupil mapping 
system that differ from the response of a conventional imaging system: an effective magni- 
fication and a distortion. These are explained in detail in Vanderbei and Traub (2005) and 
Traub and Vanderbei (2003). The magnification, in particular, is due to an overall narrowing 
of the exit pupil as compared to the entrance pupil. It is this magnification that provides 
pupil mapped systems their smaller effective inner working angle. The techniques in Section 
3 will allow us to compute the exact off-axis diffraction pattern of an apodized pupil mapped 
coronagraph and thus to see these effects. 

While the effective magnification of a pupil mapping system results in an inner working 
angle advantage of about a factor of two, it does not produce high-quaity diffraction limited 
images of off-axis sources because of the distortion inherent in the system. Guyon (2003) 
proposed the following solution to this problem. He suggested using this system merely as a 
mechanism for concentrating (on-axis) starlight in an image plane. He then proposed that 
an occulter be placed in the image plane to remove the starlight. All other light, such as the 
distorted off-axis planet light, would be allowed to pass through the image plane. On the 
back side would be a second, identical pupil mapping system (with the apodizers removed), 
that would "umap" the off-axis beam and thus remove the distortions introduced by the 
first system (except for some beam walk — see Vanderbei and Traub (2005)). A schematic 
of the full system (without the occulter) is shown in Figure 3. Note that we have spaced 
the lenses one focal length from the flat sides of the two lenses. As noted in Vanderbei and 
Traub (2005), such a spacing guarantees that these two flat surfaces form a conjugate pair 
of pupils. 



3. Diffraction Analysis 

In Vanderbei (2006), it was shown that a simple Fresnel analysis is inadequate for vali- 
dating the high-contrast imaging capabilities we seek. Hence, a more accurate approximation 
was presented. In this section, we give a similar but slightly different approximation that is 
just as effective for studying pupil mapping but is better suited to the full system we wish 
to analyze. 
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3.1. Propagation of General Wavefronts 

The goal of this section is to derive an integral that describes how to propagate a scalar 
electric field from one plane perpendicular to the direction of propagation to another parallel 
plane positioned downstream of the first. We assume that the electric field passes through 
a lens at the first plane, then propagates through free space until reaching a second lens at 
the second plane through which it passes. In order to cover the apodized pupil mapping case 
discussed in the previous section, we allow both the entrance and exit fields to be apodized. 

Suppose that the input field at the first plane is E in (x,y). Then the electric field at a 
particular point on the second plane can be well-approximated by superimposing the phase- 
shifted waves from each point across the entrance pupil (this is the well-known Huygens- 
Fresnel principle — see, e.g., Section 8.2 in Born and Wolf (1999)). If we assume that the two 
lenses are given by radial "height" functions h(r) and h(f), then we can write the exit field 
as 




E aat (x,y)=A oat (r) / / — - -e 2mQ ^ x ^ x A iQ (r)E in (x,y)dydx, (9) 

\iQ{x,y,x,y) 



— CO — CO 



where 

Q(x, y, x, y) = \J{x - x) 2 + (y - y) 2 + (h{r) - h{r)) 2 + n{Z - h(r) + h{f)) (10) 

is the optical path length, Z is the distance between the planar lens surfaces, A in (r) denotes 
the input amplitude apodization at radius r, A out (f) denotes the output amplitude apodiza- 
tion at radius f, and where, of course, we have used r and f as shorthands for the radii in 
the entrance and exit planes, respectively. 

As before, it is convenient to work in polar coordinates: 

OO 2-7T 

E out (r,9)=A out (r) [ [ ~ e 2 ^ r ' d -^A m (r)E in (r : 9)rd9dr, (11) 

J J \iQ(r,r,9 — 9) 
oo 

where 

Q(r,r,9) = \jr 2 - 2rrcos9 + r 2 + (h(r) - h(r)) 2 + n(Z - h{r) + h{f)). (12) 

For numerical tractability, it is essential to make approximations so that the integral over 9 
can be carried out analytically, thereby reducing the double integral to a single one. To this 
end, we need to make an appropriate approximation to the square root term: 



sjr 2 - 2rf cos# + r 2 + {h{r) - h{r)) 2 . (13) 
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A simple crude approximation is adequate for the 1/Q(f,r,9 — 9) amplitude-reduction 
factor in Eq. (11). We approximate this factor by the constant 1/Z. 

The Q(f,r,9 — 6) appearing in the exponential must, on the other hand, be treated 
with care. The classical Fresnel approximation is to replace S by the first two terms in a 
Taylor series expansion of the square root function about (h(r) — h(f)) 2 . As we already 
mentioned, this approximation is too crude. It is critically important that the integrand be 
exactly correct when the pair (r, f) correspond to rays of ray optics. Here is a method that 
does this. First, we add and subtract S(f, r, 0) from Q(f,r,9) in Eq. (11) to get 

Q(f,r,9-9) = S(r,r,9-9)-S(r,r,0) + S(r,r,0) + \n\(h(r)-h(r) S ) 
S(r, r,9 -9f - S(f,r, 0) 2 



+ S(f, r, 0) + \n\ (h(f) - h(r)j 



S(f,r,e-e) + S(f,r,0) 

= rf-rrcos(0-fl) /- _ \ 

(S(f,r,9-9) + S(f,r,0))/2 K J M V v; y, J K J 

So far, these calculations are exact. The only approximation we now make is to replace 
S(r, r,9 — 9) in the denominator of Eq. (14) with S(r, r, 0) so that the denominator becomes 
just S(f, r, 0). Putting this all together, we get a new approximation, which we refer to as 
the S-Huygens approximation: 

oo 2tt 

E out (r,9) « J K(r,f) j 'e 2 ^^)^^^)^^^^ (15) 



where 

K fa f ) = e 2ni( w ^+S(f,r,0)+\n\(h{f)-h(r)))/X ^ 

(note that we have dropped an exp(2ninZ / \) factor since this factor is just a constant unit 
complex number which would disappear anyway at the end when we compute intensities). 

The only reason for making approximations to the Huygens-Fresnel integral (9) is to 
simplify the dependence on 9 so that the integral over this variable can be carried out 
analytically. For example, if we now assume that the input field E in (r,9) does not depend 
on 9, then the inner integral can be evaluated explicitly and we get 

00 

E out (r,9) « 2nA ™f ] jK(r,f)J ( A / ( ~^ 0) ) 3n(0 A n {r)rdr, (17) 



Removing the dependency on 9 greatly simplifies computations because we only need to 
compute a ID integral instead of 2D. In the next subsection we will show how to achieve 
similar reductions in cases where the dependence of E in on 9 takes a simple form. 
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Figure 4 shows plots characterizing the performance of an apodized pupil mapping 
system analyzed using the techniques described in this section. The specifications for this 
system are as follows. The designed-for wavelength is 632. 8nm. The optical elements are 
assumed to be mirrors separated by 0.375m. The system is an on-axis system and we 
therefore make the non-physical assumption that the mirrors don't obstruct the beam. That 
is, the mirrors are invisible except when they are needed. The mirrors take as input a 
0.025m on-axis beam and produce a 0.025m pupil-remapped exit beam. The second mirror 
is oversized by a factor of two; that is, its diameter is 0.050m. The postapodizer ensures 
that only the central half contributes to the exit beam. The first mirror is also oversized 
appropriately as shown in the upper-right subplot of Figure 4. After the second mirror, the 
exit beam is brought to a focus. The focal length is 2.5m. The lower-right subplot in Figure 4 
shows the ideal PSF (in black) together with the achieved PSF at three wavelengths: at 70% 
(green), 100% (blue), and 130% (red) of the design wavelength. At the design wavelength, the 
achieved PSF matches the ideal PSF almost exactly. Note that there is minor degradation 
at the other two wavelengths mostly at low spatial frequencies. 

We end this section by pointing out that the S-Huygens approximation given by (15) is 
the basis for all subsequent analysis in this paper. It can be used to compute the propagation 
between every pair of consecutive components in apodized pupil mapping and concentric 
ring systems. It should be noted that the approximation does not reduce to the standard 
Fresnel or Fourier approximations even when considering such simple scenarios as free-space 
propagation of a plane wave or propagation from a pupil plane to an image plane. Even for 
these elementary situations, the S-Huygens approximation is superior to the usual textbook 
approximations. 



3.2. Propagation of Azimuthal Harmonics 

In this section, we assume that E{r, 9) = E(r)e me for some integer n. We refer to such 
a field as an nth-order azimuthal harmonic. We will show that an nth-order azimuthal har- 
monic will remain an nth-order azimuthal harmonic after propagating from the input plane 
to the output plane described in the previous section. Only the radial component E(r) 
changes, which enables the reduction of the computation from 2D to ID. Arbitrary fields 
can also be propagated, by decomposing them into azimuthal harmonics and propagating 
each azimuthal harmonic separately. Computation is thus greatly simplified even for arbi- 
trary fields, especially for the case of fields which can be described by only a few azimuthal 
harmonics to a high precision, such as Zernike aberrations, which we consider in subsection 
3.3. This improvement in computation efficiency is important, because a full 2D diffraction 
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simulation of an apodized pupil mapping system with the precision of greater than 10 10 
typically overwhelms the memory of a mainstream computer. By reducing the computation 
from 2D to ID, however, the entire apodized pupil mapping system can be simulated with 
negligible memory requirements and takes only minutes. 

Theorem 1 Suppose that the input field in an optical system described by (15) is an nth- 
order azimuthal harmonic E [n (r, 9) = E i (r)e m6 for some integer n. Then the output field is 
also an nth-order azimuthal harmonic E out (r,9) = E {f)e m6 with radial part given by 

oo 

E (r) = J K(r,r)Ek(r)J n (^ Ag( - r>0) j A- m (r)rdr. 



Proof. We start by substituting the azimuthal harmonic form of E in into (15) and regroup- 
ing factors to get 

oo 2n 

E out (rJ) = J K(r,f) j eH- t] W^)/ x E i (r)e ine d9 A in (r)rdr 

o 

K(r,f)Ei(r) [ e 2 ^~ r '^^) ,X e in{e ^ ) d9 A hi {r)rdr. 



o o 

oo 2n 

A nut (f) 



MZ 







(18) 

The result then follows from an explicit integration over the 9 variable: 

oo 


□ 



3.3. Decomposition of Zernike Aberrations into Azimuthal Harmonics 

The theorem shows that the full 2D propagation of azimuthal harmonics can be com- 
puted efficiently by evaluating a ID integral. However, suppose that the input field is not an 
azimuthal harmonic, but something more familiar, such as a (/,m)-th Zernike aberration: 
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E . m ( r ,6) = e leZ ^ r/a)cos{me \ (19) 

where e is a small number. (e/2n and e/n are the peak-to-valley phase variations across the 
aperture of radius a for m = and m/0, respectively.) 

Recall that the definition of the nth-order Bessel function is 

Jn ^ = i / elXC ° S6em(>dd - ( 20 ) 
o 

From this definition we see that i k J n (x) are simply the Fourier coefficients of e txcos ( e ) , Hence, 
the Fourier series for the complex exponential is given simply by the so-called Jacobi-Anger 
expansion 

oo 

e ixcose = i k Jk{xV k6 - (21) 

k=— oo 

The Zernike aberration can be decomposed into azimuthal harmonics using the Jacobi-Anger 
expansion: 

oo 

e ieZ r (r/a)cos(m6) = % k J k {eZ™ (r j O))^™ 6 (22) 

k=— oo 

oo 

= Jo(eZr(r/a)) + £ i k J k (eZr(r/a))e tkme (23) 
k=i 

Note that 

1^)1-^(1)" 

for < i < 1. Hence, if we assume that e ~ 10~ 3 , then the fc'th term is of the order 
10~ 3fc . The field amplitude in the high-contrast region of the PSF will be dominated by 
the k — 1 term and be on the order of 10 -3 . If we drop terms of k = 3 and above, we 
are introducing an error on the order of 10 -9 in amplitude. The error in intensity will be 
dominated by a cross-product of the k — 3 and the k = 1 term, or 1CT 12 across the dark 
region. So, in this case, Zernike aberrations can be more than adequately modeled using 
just 3 azimuthal harmonic terms. For e ~ 1CT 2 , the number of terms goes up to 5 for an 
error tolerance of 10~ 12 . In practice, even this small number of terms was actually found to 
be overly conservative. 

In order to compute the full 2D response for a given Zernike aberration, we simply 
decompose it into a few azimuthal harmonics, propagate them separately, and sum the 
results at the end. This method could also be applied to any arbitrary field. 
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4. Simulations 

The entire 4-mirror apodized pupil mapping system can be modeled as the following 
sequence of 7 steps: 

1. Propagate an input wavefront from the front (flat) surface of the first pupil mapping 
lens to the back (flat) surface of the second pupil mapping lens as described in Section 
2.3. 

2. Propagate forward a distance /. 

3. Propagate through a positive lens with focal length / to a focal plane / units down- 
stream. 

4. Multiply by star occulter. 

5. Propagate through free-space a distance / then through a positive lens to recollimate 
the beam. 

6. Propagate forward a distance /. 

7. Propagate backwards through a pupil mapping system having the same parameters as 
the first one. 

A similar analysis can be carried out for a concentric ring shaped pupil system, or even 
a pure apodization system, as follows: 

1. Choose A [n to represent either the concentric ring binary mask or some other az- 
imuthally symmetric apodization. 

2. Choose h as appropriate for a focusing lens and let h = 0. 

3. Propagate through this system a distance / to the image plane. 

4. Multiply by star occulter. 

The theorem can be applied to every propagation step, so that an azimuthal harmonic 
will remain an azimuthal harmonic throughout the entire system. Hence, our computation 
strategy was to decompose the input field into azimuthal harmonics, propagate each one 
separately through the entire system by repeated applications of the theorem, and sum 
them at the very end. 
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Figure 5 shows a cross section plot of the PSF as it appears at first focus and second 
focus in our apodized pupil mapping system (the first focus plot is indistinguishable from the 
case of ideal apodization or concentric ring shaped pupils). There are two plots for second 
focus: one with the occulter in place and one without it. Note that without the occulter, the 
PSF matches almost perfectly the usual Airy pattern. With the occulter, the on-axis light 
is suppressed by ten orders of magnitude. 

The electric field for a planet is just a slightly tilted and much fainter field than the field 
associated with the star. Hence, the methods presented here (specifically using the (1, 1)- 
Zernike, i.e. tilt) can be used to generate planet images. Some such scenarios are shown in 
Figure 6. The first row shows how an off-axis source, i.e. planet, looks at the first focus. 
As discussed earlier, at this focal plane off-axis sources do not form good images. This is 
clearly evident in this figure. The second row shows the planet as it appears at the second 
image plane, which is downstream from the reversed pupil mapping system. In this case, the 
off-axis source is mostly restored and the images begin to look like standard Airy patterns as 
the angle increases from about 2X/D outward. Figure 7 shows corresponding cross sectional 
plots for the apodized pupil mapping system at the second focus. The third row in Figure 6 
shows how a planet would appear at a focal plane of a concentric ring shaped pupil system. 

Figure 8 shows how the off-axis source is attenuated as a function of the angle from 
optical axis, for the case of our apodized pupil mapping system (at second focus) and the 
concentric ring coronagraph. For the case of apodized pupil mapping, the 50% point occurs 
at about 2.5X/D. 

Figures 9 and 10 show the distortions/leakage from an on-axis source in the presence 
of various Zernike aberrations, for apodized pupil mapping and the concentric ring shaped 
pupil systems, respectively. The Zernike aberrations are assumed to be l/100th wave rms. 

Figure 11 shows the corresponding cross-section sensitivity plots for both the apodized 
pupil mapping system and the concentric ring shaped pupil system. From this plot it is 
easy to see both the tighter inner working angle of apodized pupil mapping systems as well 
as their increased sensitivity to wavefront errors. Finally, Figure 12 demonstrates contrast 
degradation measured at three angles, 2, 4, and 8X/D, as a function of severity of the Zernike 
wavefront error. The rms error is expressed in waves. 

5. Conclusions 

We have presented an efficient method for calculating the diffraction of aberrations 
through optical systems such as apodized pupil mapping and shaped pupil coronagraphs. We 
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presented an example for both systems and computed their off-axis responses and aberration 
sensitivities. Figures 11 and 12 show that our particular apodized pupil mapping system is 
more sensitive to low order aberrations than the concentric ring masks. That is, contrast 
and IWA degrade more rapidly with increasing rms level of the aberrations. Thus, for a 
particular telescope, our pupil mapping system will achieve better throughput and inner 
working angle, but suffer greater aberration sensitivity. 

We note that there is a spectrum of apodized pupil mapping systems, out of which we 
selected but one example. The two extremes, pure apodization and pure pupil mapping, 
both have serious drawbacks. On the one end, pure apodization loses almost an order of 
magnitude in throughput and suffers from an unpleasantly large IWA. At the other extreme, 
pure pupil mapping fails to achieve the required high contrast due to diffraction effects. 
There are several points along this spectrum that are superior to the end points. We have 
focused on just one such point, which is similar to the design suggested by Guyon et al. 
(2005). We leave it to future work to determine if this is the best design point. For example, 
clearly one can improve the aberration sensitivity by relaxing the inner working angle and 
throughput requirements. Such analysis is beyond the scope of this paper, but we have 
provided here the tools to analyze the sensitivity of these kinds of designs. 
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Fig. 1. — Pupil mapping via a pair of properly figured lenses. Light travels from top to 
bottom. 




-0.5 0.5 -60 -40 -20 20 40 60 



Pupil-Plane Radius in fraction of Aperture Image-Plane Radius in L/D radians 



Fig. 2. — Left. An amplitude profile providing contrast of 10 10 at tight inner working 
angles. Right. The corresponding on-axis point spread function. 
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Fig. 3. — The full pupil mapping system includes: a pair of lenses to shape the amplitude 
into a prolate-spheroidal-like amplitude profile, a focusing lens that concentrates the on- 
axis starlight into a small central lobe where a (not-depicted) occulter can block this light, 
followed by a recollimating lens and finally a reverse pupil mapping system that reforms the 
pupil with the starlight removed but any planet light (if present) intact. This final pupil is 
then fed a final focusing element (not shown) to form an image of off-axis sources. 
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2nd Pupil Amplitude Map 
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Fig. 4. — Analysis of an apodized pupil mapping system using the S-Huygens approxima- 
tion with z = 15.D and n = 1.5. Upper-left plot shows in red the target high-contrast 
amplitude profile and in blue the amplitude profile computed using the S-Huygens approx- 
imation through the apodized pupil mapping system. The other two curves depict the pre- 
and post-apodizers. Upper-right plot shows the lens profiles, red for the first lens and blue 
for the second. The lens profiles h and h were computed using a 5, 000 point discretization. 
Lower-left plot shows in red the computed optical path length Qo(f) and in blue the phase 
map computed using the S-Huygens propagation computed using a 5, 000 point discretiza- 
tion. Lower-right plot shows the PSF computed at three different wavelengths; the design 
value, 30% above that value, and 30% below it. 
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Fig. 5. — On-axis PSF at first focus (before occulter) and at second focus for cases of with and 
without occulter. Without the occulter, the second-focus PSF almost perfectly matches the 
usual Airy pattern. However, with the occulter, the second-focus on-axis PSF is suppressed 
by ten orders of magnitude. 
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Fig. 6. — 2D pictures of planets for apodized pupil mapping and concentric rings. First 
row shows 2D intensity plots at first focus behind the occulter for planets at various angles 
relative to the on-axis star. Note that the system fails to form a clean image of the planets. 
Second row shows analogous plots at second focus. Note that the wavefront for the off-axis 
planet is mostly restored and the images begin to look like standard Airy patterns as the 
angle increases. Third row shows analogous plots for a concentric ring mask. 
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Fig. 7. — Cross-sectional plots from the second row plots in Figure 6. Note that for angles 
of 3A/.D and above, the restored PSF looks very much like an Airy pattern with very little 
energy attenuation. However, as the angle decreases, the pattern begins to distort and the 
throughput begins to diminish. 
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Fig. 8. — Off-axis source attenuation as a function of angle. Note that the 50% point occurs 
at about 2.5X/D. 
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Fig. 9. — Apodized pupil mapping sensitivities to the first nine Zernike aberrations. In 
each case, the rms error is l/100th wave. The plots correspond to: Piston (0,0), Tilt (1, 1), 
Defocus (2,0), Astigmatism (2,2) Coma (3,1), Trefoil (3,3), Spherical Aberration (4,0), 
Astigmatism 2nd Order (4,2), Tetrafoil (4,4). 
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Fig. 10. — Concentric-ring mask sensitivities to the first nine Zernike aberrations. 
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Fig. 11. — Radial profiles associated with the previous two Figures and overlayed one on the 
other. The green plots are for apodized pupil mapping whereas the blue plots are for the 
concentric ring mask. 
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Fig. 12. — Contrast degradation measured at three angles, 2, 4, and 8X/D as a function of 
severity of the Zernike wavefront error measured in waves. 



